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Entanglement distribution between trapped-atom quantum memories, viz. single atoms in optical 
cavities, is addressed. In most scenarios, the rate of entanglement distribution depends on the 
efficiency with which the state of traveling single photons can be transferred to trapped atoms. This 
loading efficiency is analytically studied for two-level, V-level, A-level, and double-A-level atomic 
configurations by means of a system-reservoir approach. An off-resonant non-adiabatic approach 
to loading A-level trapped-atom memories is proposed, and the ensuing trade-offs between the 
atom-light coupling rate and input photon bandwidth for achieving a high loading probability are 
identified. The non-adiabatic approach allows a broad class of optical sources to be used, and 
in some cases it provides a higher system throughput than what can be achieved by adiabatic 
loading mechanisms. The analysis is extended to the case of two double-A trapped-atom memories 
illuminated by a polarization-entangled biphoton. 

PACS numbers; 03.67.Hk, 03.67.Mn, 42.50.Ct, 42.50.Pq 
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I. INTRODUCTION 



The future of quantum information science — which 
holds the promise of revolutionary improvements in com- 



' prom: 

putation [ll, [2|, secure communication Q , and precision 
measurement [J| — depends on our ability to develop its 
demanding technology. Some of the required building 
blocks have already been developed, e.g., first genera- 
tion quantum key distribution systems are now commer- 
cially available Q . With recent advances in single-photon 
technology — the development of high-flux sources of en- 
tangled photons @, 01 as well as efficient (photon- number 
resolving) single-photon detectors [1, [1,1131 — is foresee- 
able that the next generation of secure long-distance com- 
munication systems will be within reach in a few years. 
The key element, which is still missing from such systems, 
is a reliable quantum memory that can store and manipu- 
late quantum information. Whereas trapped-ion systems 
are promising choices for quantum processing, the lack of 
an efficient optical interface hinders their being used in 
long-distance quantum communication systems. Some 
neutral alkali atoms, however, e.g. rubidium and cesium, 
are appropriate choices for capturing the states of single 
photons and storing them in their metastable hyperfine 
levels. In order to enhance the interaction of single pho- 
tons and atomic systems, we have to either trap single 
atoms in high-finesse optical cavities fill, [l^ Ha . [T^ or 
use ensembles of atoms [H [H, [13, [lllia llj In ei- 
ther case, it is an interesting problem to quantify how 
efficiently we can transfer the state of a single photon 
to an atomic quantum memory, or, in other words, how 
efficiently we can load a quantum memory with a de- 
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sired state. In this paper, we study the loading problem 
in the context of trapped-atom quantum memories, viz. 
single atoms trapped in high-Q optical cavities, and find 
analytical solutions for the loading probability in differ- 
ent scenarios in which either a single photon or a pair 
of entangled photons are driving the quantum memory 
units. That analysis results in a novel non-adiabatic ap- 
proach to loading quantum memories with A and double- 
A atomic configurations. We also investigate the effect 
of various system parameters on the loading probability, 
including the atom-light coupling rate, the input light 
bandwidth, the cavity/atomic decay rate, and timing off- 
sets. Atomic-ensemble memories will not be explicitly 
addressed in this paper, however, most of the analytical 
results that we obtain are directly applicable to such sys- 
tems. We elaborate more on this issue in the discussion 
section. 

Entanglement is a quantum resource by which two 
parties share a joint state that cannot be written in 
a tensor-product form. This state provides a stronger- 
than-classical correlation between two quantum systems 
by which one can perfectly infer quantum measurement 
results made on one system by observing the results of 
quantum measurements made on the other. This quan- 
tum correlation is behind various applications in quan- 
tum communication [2l| and distributed quantum pro- 
cessing [l^l. In teleportation, entanglement serves as a 
quantum wire. Once established, it can be used for one- 
time communication between the transmitter and the re- 
ceiver. In distributed quantum processors, entanglement 
serves as a computational resource. The main trick in 
both of these applications is that we have connected our 
parties, in advance of any subsequent actions, via entan- 
glement. This way, we need not worry about the channel 
loss or physical transportation of a qubit once we have 
established entanglement. We ought to worry, however, 
about how to generate, distribute, and maintain entan- 
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FIG. 1: (Color online) Three architectures for entanglement 
distribution: (a) To-the-memory architecture, in which two 
entangled photons, generated at the source, carry the entan- 
gled state to the quantum memory units and load them with 
this state, (b) From-the-memory architecture, in which, by 
the use of a common source, each memory generates a fly- 
ing qubit entangled with itself. Entanglement swapping is 
then accomplished by means of a Bell-state measurement on 
the photons, at the midpoint, which leaves the memories in 
an entangled state, (c) Memory-to-memory conflguration, in 
which a flying qubit entangled with one of the memories is 
generated that then propagates to and loads the other mem- 
ory, i.e., transfers its state to that memory. 



glement over long distances, i.e., distances over which our 
physical system undergoes decoherence. 

There are three commonly suggested approaches for 
long-distance entanglement distribution, all of which use 
photons for the transmission of quantum information and 
atomic systems for its storage. They are shown schemat- 
ically in Fig. [TJ In the first approach, we produce a pair 
of entangled photons at an optical source, and let the 
photons travel to and be stored in the quantum mem- 
ory units [23'! . In the second approach, we use entan- 
glement swapping (23 | to distribute the entanglement. 
We first generate entanglement between a standing qubit 
and a flying qubit at each end of Fig. [Ijb). We then let 
the photons propagate to the midpoint between memo- 
ries where a Bell-state measurement apparatus destroys 
them, lea ving the quantum memory units in an entangled 
state [TH. l25l| The last architecture, shown in Fig. [Ijc), 
is a hybrid of the previous ones. In this approach, we 
generate a single photon entangled with one of the quan- 
tum memories, and let this photon propagate to and load 
the other memory, i.e., transfer its state to that memory 

In all of these entanglement distribution schemes, a 
figure of merit for system performance is its throughput, 
viz. how often we obtain an entangled pair of quantum 
memories. Throughput is a function of different system 
parameters. For instance, in the to-the-memory archi- 
tecture, Fig. [Ija), it depends on the source generation 
rate of entangled photons, the path and coupling losses, 
and the loading efficiency. System loss has the unavoid- 
able effect of reducing the throughput. Loss, however, 
does not affect the source generation rate or the nature 
of interaction between a single photon and a single atom. 



Moreover, for trapped-atom quantum memories, there 
are nondestructive loading verification techniques that 
allow us to detect photon loss Hence, throughout 

the paper, we assume that the system is lossless, and we 
focus on the interaction between the source output and 
the memory modules. 

The interaction between a trapped atom and light is 
governed by several factors, among them are the input 
photon bandwidth, the cavity decay rate, and the atom- 
light coupling rate. In order to realize a strong atom-light 
interaction, the photons must have bandwidths compara- 
ble to atomic linewidths, which are typically on the order 
of tens of MHz. Higher optical bandwidths may afford 
increased source generation rates — as is the case for spon- 
taneous parametric downconversion sources whose band- 
widths are on the order of THz [13, S [H, H^]— and, 
hence, increased total throughput. However, if higher 
source bandwidth comes with lower spectral brightness, 
throughput may be decreased, owing to a decrease in 
the loading probability. It may be possible to main- 
tain a high loading probability by employing a higher 
atom-light coupling rate. This last approach requires a 
more demanding implementation, however, because, for 
trapped-atom quantum memories, the coupling rate is 
increased by reducing the cavity length. All these is- 
sues necessitate developing a quantitative loading analy- 
sis that quantifies the interplay between the input pho- 
tons' bandwidths and the quantum memory parameters, 
i.e., cavity decay rate and atom-light coupling rate. This 
paper seeks to answer these questions, which are of prac- 
tical importance in designing quantum memory systems, 
by quantifying the trade-offs inherent in adiabatic and 
non-adiabatic loading mechanisms. 

Another important issue in dealing with the interac- 
tion of trapped atoms and light is decoherence. There 
are a variety of mechanisms by which our system may 
decohere. Decay of photons out of the cavity and atomic 
spontaneous emission are the two factors that will be con- 
sidered in our analysis. The latter effect can be alleviated 
by using off-resonant transitions, which will be employed 
in many scenarios in this paper. There are other sources 
of decoherence — some related to the employed trapping 
and cooling schemes [Slj , and some due to the oversimpli- 
fied picture that we use for atomic systems [33 | — which 
will not be addressed in our paper. 

The rest of this paper is organized as follows. In Sec- 
tion |TT1 we study the loading process for single trapped- 
atom quantum memories that are illuminated by single 
photons. The methodology that we employ in this sec- 
tion will be applicable to systems that use the memory- 
to-memory configuration for entanglement distribution. 
We start with the basic case of a two-level atom, and then 
extend our results to A- and double-A-level atoms. We 
propose a non-adiabatic loading mechanism for the latter 
cases, and compare it to the previously proposed adia- 
batic loading techniques [1^, |33| . In Section IIIII we treat 
an example of the to-the-memory configuration in which 
the source output is a polarization-entangled biphoton. 
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FIG. 2; (Color online) A single two-level atom trapped in 
a single-ended high-finesse optical cavity. A single photon 
illuminates the cavity at the proper frequency to excite the 
atom to its upper state. 



and the trapped atoms have double-A configurations. 
This is a close approximation to the architecture pro- 
posed by researchers from the Massachusetts Institute of 
Technology and Northwestern university (termed MIT- 
NU hereafter) [ll|, [H, [111. Thus the results presented 
here comprise the first loading analysis of the MIT-NU 
system that accounts for the intracavity atoms. 



II. TRAPPED-ATOM QUANTUM MEMORIES 
DRIVEN BY SINGLE PHOTONS 



In this section, we evaluate the loading probability for 
trapped-atom quantum memories, viz. single atoms in 
high-finesse optical cavities, driven by single photons. 
The loading protocol is successful if we transfer the pho- 
ton state to the corresponding atomic levels. We do 
not deal with the complexity associated with trapping 
neutral atoms, and assume that the single atoms can be 
trapped at a fixed point in the cavity for a time sufficient 
for loading to occur [l2[. Throughout the paper, opti- 
cal cavities are assumed to be lossless and single ended. 
The external photon illuminates the cavity through its 
partially reflecting mirror in a specific spatial and polar- 
ization mode matched to the cavity's mode of interest. 
We assume that the incoming light has such a narrow 
bandwidth that it only interacts with a single longitudi- 
nal mode of the cavity. To model the temporal content 
of the incoming photon, we consider a reservoir of har- 
monic oscillators corresponding to a continuum of modes, 
as suggested in [35] . This is also in accord with the quan- 
tum representation of traveling- wave light |36| . We then 
employ a system-reservoir approach to analyze the load- 
ing dynamics. 



A. Two-level Trapped Atoms 

Suppose that there is a single two-level atom, with 
excited state |e) and ground state \g), inside a single- 
ended high-Q cavity with decay rate k. We assume that 
the frequency-Wa atomic transition between |e) and \g) is 
coupled, at coupling rate g, to the cavity field operator 
b of frequency ujq . This implies a detuning A = loq — uja, 
which is assumed to be much less than loq. In general, g 



and K can both be functions of time. Here, we assume 
that K is a constant. Suppose that this trapped-atom 
module, with the atom initially in its ground state and 
no photon in the cavity, is illuminated by a single photon, 
as shown in Fig. [2l This driving source can be modeled 
by a set of annihilation operators , each corresponding 
to a different temporal (spectral) mode of frequency w, 
that satisfy [a[^,a^,] — S{uj — uj') We assume that 

the driving field is initially in the following single-photon 
state 



(1) 



where Jdtj|0(a;)p = 1, and jl^^) = aj^|0)fl is the 
multi-mode state representing one photon at frequency 
Lo and the vacuum state for all other modes. Here, 
\0) R is the reservoir's multi-mode vacuum state. We 
can equivalently describe the initial state in the time do- 
main by introducing the temporal pulse shape ^(t) — 
J du!e^'^'^*<j>{uj)/\/2n associated with the incoming pho- 
ton, and single-photon states |lt) = / daje'"'|lt^)/\/27r 
associated with time t, as shown in Eq. ([T]). Here and 
elsewhere in this paper, we assume that the initial time 
is before any possible interaction between the source and 
the memory. For simplicity, and without loss of general- 
ity, we assume that this initial time is 0. We will relax 
this assumption after obtaining our final results. 

The Hamiltonian for the above system is given by [s^l 

Ha = Hcc + fuuJad-ee + %(&'^(Tge + &(7eg), (2) 

where ct^ = \i) {j\, i,j e {g,e}, and 
Hcc = fi / dijj udljduj + huJoo'b + hr / duj{aljb + b'' di^), 



(3) 

with r = k/tt being the coupling constant that con- 
nects the external world to the cavity [sH]. Because of 
the source's narrow bandwidth, we can and will assume 
that all integrals in Eq. ([3|) run from — oo to +00. 

Neglecting for now the decoherence mechanisms that 
influence the time evolution, the quantum system consist- 
ing of the cavity, the atom, and the reservoir is closed, 
with exactly one excitation at any time t. A general 
quantum state of the system at time t can thus be writ- 
ten as 

\m) = J dcjaUt)\W\0)b\9)+e-''^"'p{tmR\l)b\g) 
+e-*-°*Ce(t)|O)fl|0)b|e), (4) 

where |fc){, represents the fc-photon Fock state of the cav- 
ity mode. 

The goal of this section is to provide analytical results 
for the memory loading probability — the probability of 
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absorbing the external photon by the trapped atom — 
|ce(t)p. The atom in our system is connected to the ex- 
ternal single photon via the cavity mode. The interaction 
of the reservoir and the cavity mode can be analyzed by 
applying the Schrodinger equation, ih\ip(t)) = Ha\^p{t)), 
to Eq. to obtain 

a^t) = -i{coauj{t) + re-''^°*/3(t)), (5a) 
j3{t) = -iEe^^o* / duj a^{t). (5b) 

Solving for a^{t) in Eq. ([5a|) . and then plugging the result 
into Eq. (j5bp . we obtain 



(6) 



P(t) = -igce{t) ~ iV2K<Pb{t) - Ki3{t), 



where $&(i) = e*'^''*$(i) is the baseband input pulse 
shape provided that the center frequency of the incom- 
ing light is wq. The second term on the right-hand side 
of the above equation accounts for the effect of the in- 
put photon. The last term accounts for the decay of the 
photon out of the cavity, and in essence, is similar to the 
decay term in the Weisskopf-Wigner theory of sponta- 
neous emission [37|, in which a two- level system initially 
in its excited state interacts with a reservoir initially in 
its vacuum state. In this way, we can incorporate the 
overall spontaneous decay, with rate 7, from the excited 
state of the atom to the modes other than the cavity's 
mode of interest, and thus obtain 



Ce{t) = iACe{t) - igl3{t) - JCe{t), 



(7) 



where the first two terms on the right-hand side of Eq. ([7]) 
represent the coherent evolution of Ce(t) in accord with 
the Schrodinger equation ih\ip{t)) = Ha\ijj{t)). 

Equations (O and ([7]) model the loading dynamics of 
a two-level trapped-atom memory illuminated by a sin- 
gle photon with pulse shape $(i). For the general case 
of time-dependent system parameters, these equations 
can be solved numerically. For the special case of time- 
independent g, K, 7, and A, we can solve the above 
equations analytically using the Laplace-transform tech- 
niques, which yield 

P^t) = f dT$fc(r)(Ac;e-"+(*-") - K'_e-«-(*-^)) 

? Jo 

(8) 

and 



3(0 = 



g\/2K 

~rjo 



dT$b(T)(e-''+(*-^) 



where 



K± = (k + 7' ± 0/2, k'± = K± 



-K,-{t-T) 



), (9) 



(10) 



The lower limit zero in the above integrals represents the 
initial time of interaction. It can be changed to accom- 
modate the input pulses that start before i = as well. 



Note that Eqs. (O and ([7|), with initial values Ce(0) = 
/3(0) — 0, are linear in response to the input pulse shape. 
For instance, assuming that $b(i) = <&i(i) + *2(0> /3(i) 
will be equal to fii{t) -\- (32{t), where /3i{t), for i — 1,2, 
is given by Eq. ([8]) with ^bit) replaced by $i(t). Simi- 
larly, we can think of the initial state in Eq. (IT|) as a su- 
perposition of its infinitesimal constituents da;0(a;)|l(^), 
to which, from Eq.([5]), we can associate slowly- varying 
probability amplitudes Ce,cj(i) given by 



/q dTdtj0(cj)e-'('^-'^o)^(e-'"+(*-^) - e-''-^*-^)). (11) 

We then have Ce{t) = / du)Ce^uj{t) as obtained in Eq. 
We use this result later in Section [1111 

The system-reservoir calculation gives a more compact 
form for the loading probability, |ce(i)P, than what we 
obtained in [3^ via the Heisenberg-Langevin equations 
[3^ for the case of on-resonance illumination at 7 = 0. 
The two results are nevertheless equivalent, and their 
equivalence can be verified by a tedious algebraic manip- 
ulation that we shall omit. 



1. Numerical results 

There are several physical parameters that are of prac- 
tical importance in loading a trapped-atom quantum 
memory: cavity decay rate, spontaneous decay rate, 
atom-light coupling rate, and input bandwidth. How- 
ever, from Eq. it can be seen that the loading proba- 
bility is governed by the ratios of these parameters. Thus, 
for A = 0, Eq. ^ can be written in terms of the dimen- 
sionless parameters kT, g/ k, 7/g, and t/T, as follows 

Ce{t) = -jj^ dM$'t(u)x 

-KT(K+/K)(t/T-tl) _ g-KT(K_/K)(i/T-tl) ^^2^ 

where T is the effective width of ^b{t), and ^b{t) = 
^'i^{t/T)/\/T, i.e., ^b{u) is a compressed/stretched ver- 
sion of <&b(i), which has unity width, in its normal- 
ized time coordinate, and has been normalized to satisfy 
Jdu\mu)\^ = l. 

For a constant value of g, what we expect from the 
atom-photon interaction is a damped Rabi oscillation. 
The incoming photon is initially transferred to the cav- 
ity mode, which interacts with the atom and drives a 
|g)-to-|e) transition. That will be followed by the re- 
verse transition which releases a photon into the cav- 
ity. These oscillations continue until the photon leaks 
out of the cavity or the atom spontaneously emits a pho- 
ton into a non-cavity mode. Figure [3] shows the loading 
probability for the two-level atom of Fig. [5] in the ideal 
case of 7 = A = for a hyperbolic secant pulse shape 
<I>^(u) — \/2sech[4(w— 1)]. In this figure, we have fixed the 
value of kT to 2 and varied the value of g/n. Counterintu- 
itively, it can be seen that increasing the coupling rate g 
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FIG. 3: (Color online) Loading probability versus normalized 
time for kT = 2 and a hyperbolic secant pulse shape at 7 = 
A = 0. Upon arrival of the photon, there is a peak in the 
loading probability, which occurs at the time that the photon 
is most likely to have absorbed by the atom. This maximum 
loading probability is a function of q/k, and there exists an 
optimum value for this parameter that maximizes the chance 
of loading. 
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FIG. 4: (Color online) The optimum coupling rate, and the 
maximum loading probability achieved at this rate, versus k,T 
for loading a two-level atom by a single photon with hyper- 
bolic secant pulse shape on resonance. Several spontaneous 
decay rates, 7, are considered. 



does not necessarily improve the loading probability, but 
there exists an optimum value of g/K, which maximizes 
the loading probability. Higher ratios than this optimum 
value just increase the number of Rabi flops between the 
two atomic levels, making it harder to find the atom in 
its excited state with high probability. However, by this 
method we cannot hold the atom in its excited state un- 
less we drive the value of g to zero at an appropriate point 
in the process, t = T^oad, when the loading probability 
attains its maximum. In the next section, we show how 
we can turn off the coupling in a A-level atom, in which 
we can use a control field to vary g. 

Figure |4] shows the optimum value of g/n versus kT 
for several values of 7 and A = 0. Several things are no- 
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FIG. 5: (Color online) Left: loading probability at kT = 2, 
g/n — 1, and 7 = A = for various pulse shapes. Right: the 
baseband pulse shapes. In order of their maximum loading 
probabilities they are: the hyperbolic secant (with the highest 
peak), rectangular, exponentially-rising, and exponentially- 
decaying pulse shapes. 



ticeable in this figure. First, for a fixed k, it can be seen 
that the required value for the coupling rate will decrease 
as we increase the pulse width, or equivalently, when 
our source's bandwidth decreases. The decrease that we 
achieve, however, by doing so is exponentially smaller as 
we go to larger values of kT. For instance (not shown on 



the graph), to get {g/K) 



optimum 



0.1, we need kT to be 



approximately 100. It can also be seen that the optimum 
coupling rate will not change significantly in the pres- 
ence of decoherence. Nonzero spontaneous decay rates, 
however, will reduce the maximum loading probability — 
which could be above 90% otherwise — as shown in the 
lower part of Fig. HI Longer pulses, which require longer 
loading times, are more exposed to the decay process, and 
so, their maximum loading probabilities are lower. Spon- 
taneous decay can be alleviated by using large enough 
detunings in a A-level atomic system. 

Figure [5] shows the loading probability as a function of 
time for g/n = 1 and kT = 2, for different pulse shapes of 
approximately the same effective width and 7 = A = 0. 
The four different shapes we have used are: a hyperbolic 
secant pulse; a rectangular pulse; an exponentially-rising 
pulse; and an exponentially-decaying pulse. The differ- 
ence between their loading probabilities is seen to be very 
minor. It seems that having matched the input band- 
width to the cavity parameters g and k, we can achieve 
high loading probability regardless of the photon pulse 
shape. For the rest of this paper, we only consider the 
hyperbolic secant pulse shape in our calculations. 



B. A-level Trapped Atoms 

The two-level atom that we analyzed in the previous 
section, while the easiest to study, is not a practical so- 
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FIG. 6: (Color online) A A-level trapped atom illuminated 
by a single photon. The cavity mode corresponding to the 
single photon drives the atom from its ground state \g} to 
the auxiliary state |r). A z-polarized beam then shelves the 
atom in the metastable level |e). Here, Qc is the vacuum Rabi 
frequency for the |(7)-to-|r) transition, and Q is the Rabi fre- 
quency associated with the control field, which is proportional 
to \Ez\. These transitions are off-resonant by detunings Ai 
and A2, respectively, as defined in the text. The control field 
may also include a time- varying phase 4iz{t). 



all decoherence mechanisms, with no initial excitation in 
the cavity and assuming that the atom is initially in its 
ground state l^), the most general time evolution for the 
state of the system takes the following form 



m)) 



duja^{t)\l^)\Q)b\g) 



m\G)+e 



MR) 



where 



|G) EE \0)n\l),\g), 
\R) EE |O)fi|0)b|r), 
\E) EE \O)n\0)b\e). 



(14) 



(15a) 
(15b) 
(15c) 



Applying the Schrodinger equation ihlip^t)) = Hi,\tp{t)) 
to Eq. dUl), we get 



lution for quantum storage. For long lifetime storage, 
we have to store the quantum information in metastable 
atomic states that have low energy levels as well as low 
decay rates, e.g. the hyperfine levels in the 5^5*1/2 or- 
bital of the rubidium atom. Whereas the direct transition 
between such metastable levels are usually dipole forbid- 
den, we can connect these states, via a third excited level, 
by inducing a Raman transition driven by a single pho- 
ton on one leg and a classical control field on the other. 
The spontaneous decay from the upper state can then 
be mitigated by including large detunings. The resulting 
A-level atom is at the core of all interesting neutral-atom 
quantum memory units. This section studies the loading 
problem in the context of trapped-atom memories with 
A-level configurations. 

Figure [6] shows a trapped A-level atom and its corre- 
sponding driving beams. Here, we assume that an ex- 
ternal single-photon beam, in the initial state given by 
Eq. ([1]), is spatially matched to a cavity mode with an- 
nihilation operator b and resonance frequency loq. This 
cavity mode drives the |g)-to-|r) transition with coupling 
rate gc- There may exist a detuning Ai = c^o — ujgr 
in this transition, where cjy denotes the transition fre- 
quency between levels \i) and for i,j g {g,r,e}. The 
second beam is assumed to be a z-polarized classical wave 
with frequency tOz and a possibly time-dependent phase 
<pz{t). This field's amplitude, which is under our control 
and may vary with time, determines the Rabi frequency 
il{t). It couples \r) and |e) via the following Hamiltonian, 
which is obtained under the rotating- wave approximation 



ngdb^ 



Gar 



CFrgb) 

+ /in(t)(e-*["^*+*^(*)la,e + e'["^*+*^'*^l<Ter), (13) 



where il{t) and (t>z{t) are assumed to be real. The \r)- 
to-|e) transition may also be off- resonant, by a detuning 
A2 = oJz—iiJer- Similar to the previous section, neglecting 



(3{t) = ~igcCr{t) - i\/2^<I>fc(t) - Kl3(t) 



(16a) 



c^{t) = ~igcl3{t) + iAiCr(t) - in{t)ce{t) - 7^0^(0 (16b) 

Ce(t) = -in{t)Cr{t) - l{A2 - Ai + 0^(i))Ce(t), (16c) 

where we used Eq. ^ in the first equation and included 
a non-cavity decay rate 7^ for the upper state. We as- 
sume that the corresponding decay rate to the state |e) 
is negligible for the purpose of loading; it comes into play 
if we need to determine the storage time of the quantum 
memory, which is not a topic of interest in this paper. 
For large enough detunings, i.e., Ai, A2 3> ^,gc, we can 
adiabatically eliminate the upper state \R) by assuming 
that Cr{t) — 0. As a result, we obtain 



Cr{t) 



9c/A 



-Pit) 



n{t)/A 



-Ce{t) 



(17) 



where F^ = 1 -f ijr/Ai. Plugging the above equation 
into Eq. (HH), we then obtain 



= /Ai..,. .git) 

-Pit) -^-^Ce(t) 



At), 



1 r 



-i{A2- Ai+^zit))Ceit), 



(18a) 



(18b) 



where g{t) — g^Slit) / Ai is the effective coupling rate be- 
tween \g) and |e) in the absence of decay. These equations 
resemble Eqs. ^ and ([7]), which govern the evolution of 
a two-level system. To see this equivalence, we define 
new variables /3'(i) = exp[j5^t/(AiFr)]/3(t) and d^it) ee 
exp[i(7^</(AiFr)]ce(i). Assuming that 7r-/Ai <C 1, and 
thus approximating |Frp by one, we then obtain 

P'[t) = -i!Mc'^{t) - iV2^$^t) - K(3'{t), (19a) 
m = -*#/3'(0 + *A4(t) - 74(0, (19b) 



where 



shape <^b{t) = v/2/Tscch(4</r), it turns out that [SE 



7 = 9{thr[^{t)/gc 



Ai- A2-0.(i), 
-5c/^^(0]/Ai. 



(20a) 
(20b) 
(20c) 



Equations (|19ap and (|19bp are similar to Eqs. ([5]) and ([7]) 
if we use the new values for 7 and A as given by Eq. ((20)) . 
and replace g and $b(i), respectively, with g{t)/Tr and 
$'{,(t). The loading probability is then given by |ce(i)p w 



e -i-c "■v"i|c^(t)|2^ Note that the term e*fc*/Ai $[,(t) 
induces a constant frequency shift to the incoming light. 
We can tilt the center frequency of the input photon from 
its original value luq to compensate for this induced shift. 
From now on, we assume that this compensation has been 
employed. 

For our loading problem, we are looking for particular 
pulse shapes for n{t) and (j)z{t) that help the \g)-to-\e) 
transition occur with high probability. The traditional 
solution to this problem employs adiabatic-transfer tech- 
niques p6l . [33L [39L [ioj , in which we slowly change the con- 
trol field so that the system follows its dark-state eigen- 
state. In our case, for Ai = A2 + 0z(t) and 'jr = 0, the 
system's dark state is given by, 



\D) 



-cosd{t)\G) +smeit)\E). 



(21) 



where cos6'(t) = f2(t)/f2o, smd{t) = gc/^o, and flo ^ 
^^■^{t) + gl- This state has the desired property that 
for n{t) > <?e, \D) « -|G), but for n{t) < g^, \D) « 
\E). So, if we start with a high value of Vl when the 
photon arrives, and then slowly reduce to zero, we can 
adiabatically transfer the system from \G) to \E). The 
timing is of crucial importance, because if we turn off the 
pump either before or long after the photon arrives we 
lose the chance of absorbing the photon. There is also a 
slight chance of jumping into states that are orthogonal 
to i.e., \R) or \B) = sm.e[t)\G) + cos 9 {t)\E) , in which 
case, the loading process has completely failed. 

Under the dark-state condition, in which the probabil- 
ity of jumping to either \R) or \B) is near zero, we can 
show that our loading problem is equivalent to loading 
a cold cavity, viz. a cavity with no atom inside, with a 
time-varying decay rate kcos^ 9{t). This can be shown 
by paralleling the derivation used in [11] and shall be 
omitted. To maximize the loading probability, we need 
to find an optimum assignment for cos9{t) that maxi- 
mizes the loading probability and satisfies our adiabatic 
conditions. In order to achieve a maximum transfer of 
free-field photons into the cavity mode, we need to min- 
imize the outgoing field components by destructively in- 
terfering the directly refiected and the circulating fields. 
Fleischhauer et al. have shown that a necessary condition 
for destructive interference is (33l |: 

- — In cos^^(<) + — ln<I'b(t) — Kcos 9{t). (22) 

{A/ L \JLL 

For our particular example of a hyperbolic secant pulse 



5cSech(4t/T) 



^[1 + tanh(4t/r)][tanh(4t/r) + kT/2 - 1] 



(23) 

In order for this Vl{t) expression to yield a positive real 
Rabi frequency, we must have kT > 4. This, on the other 
hand, implies that for a successful adiabatic transfer, our 
input pulse must be long enough so that we can slowly 
change the quantum memory state. 

Here, we propose a simple non-adiabatic approach for 
loading a A-level atom, which does not impose any re- 
strictions on the input pulse shape and does not need any 
adiabatic pulse shaping for the control field. Our method 
is based on what we observed in Section III Al for a two- 
level atom with a constant coupling rate g. There, we 
realized that a maximum loading probability of greater 
than 90% was achievable, provided that we could turn 
off the atom-light coupling at t = TLoad- Here, we show 
that this is indeed possible to do for a A-level atom, 
by applying a constant control field for an appropriate 
finite-duration time interval. In other words, we assume 
(j)z{t) = and take fl{t) to be 



n{t) = 



0, 



t < Tlo 
t>TLo 



(24) 



From Eq. (fT9|) . the effective coupling rate, g = gci^/Ai, is 
proportional to the control field's amplitude. By turning 
off the control field at i = TLoad, the couphng rate g 
vanishes for t > TLoad- Hence, if we are in the state |e) 
at t = TLoad, we will stay there until a decay process 
returns the atom to its ground state. Solving for c'^{t) in 
Eq. ([I9b| using Eq. ([9]), the loading probability is then 
given by 



|Ce(t)P 



2Kg2 

w 



dre-^'''^'-(*-^)/'^i$6(T)(e-'=+(*~^) - e-'"-^*-^)) 



for t < TLoad, 



(25) 



where k± and ^ are given by Eq. (jlOp using 7 and A 
as given by Eq. ([20|l . The loading probability is then 
given by \ce{TLoad)\'^ - The effective detuning, A = (g^ — 
ri^)/Ai -I- Al — A2, can be forced to zero by properly 
choosing Ai and A2, thus compensating for the induced 
Stark shift and enabling us to obtain a higher loading 
probability. 



1. Adiabatic versus non-adiabatic loading: numerical 
comparison 

In this section, we compare our proposed non-adiabatic 
scheme to loading mechanisms that use adiabatic- 
transfer techniques. Adiabatic transfer is a well-studied 
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problem in the literature. To make our comparison ex- 
plicit, we will consider the particular examples of adia- 
batic loading that have been proposed in [H] and [33l |. 
In [26*1, the authors have devised a method for trans- 
ferring the state of a trapped-atom quantum memory 
to another trapped-atom quantum memory. Their ap- 
proach is based on adiabatically transferring the state of 
one memory to a single photon, which then propagates 
to and loads the other memory. Their loading process is 
facilitated by forcing the photon's pulse shape to be sym- 
metric, so that the receiver can employ a time-reversed 
version of the control pulse that was used at the trans- 
mitter. The desired pulse shape for the control field, un- 
der the dark-state condition at zero effective detuning, 
can be found numerically by solving the corresponding 
Schrodinger equations. However, the approach from |26j| 
is not suitable for an incoming photon with an arbitrary 
pulse shape. In (s^, the authors have employed the adi- 
abatic transfer to load an atomic ensemble with the state 
of a single photon. They have employed the dark-state 
approximation under the two-photon resonance condi- 
tion with a constant-phase control field. Their approach 
sets certain limitations on the length of the input pulse 
shape. In fact, in order to fulfill the dark-state condition, 
the input pulse shape must be longer than a threshold 
value. Unlike these two adiabatic mechanisms, our non- 
adiabatic approach puts no constraints on the input pulse 
shape. 

Several issues make our work in this paper, and the 
numerical results presented in this section, distinct from 
previous work on loading quantum states into neutral 
atoms. The first issue is our accounting for a nonzero 
probability for populating the bright state \B). We will 
see in this section how the nonzero probability of being 
in the bright state affects the loading performance. That 
also makes it possible to scrutinize the dependence of the 
loading probability, in both adiabatic and non-adiabatic 
mechanisms, on the key system parameters, e.g., gc/i^ 
and kT . This evaluation is one of the original contribu- 
tions of this paper. Finally, our use of a constant control 
field allows us to obtain analytical results for loading dy- 
namics in the non-adiabatic case. 

It is of practical importance to know at what values 
of input bandwidth, represented by kT, and atomic cou- 
pling rate, represented by gc/i^, a desired system perfor- 
mance can be achieved. For our non-adiabatic approach, 
at 7r = A = 0, system performance is governed by the ef- 
fective coupling rate g = gcfl/Ai. The optimum value of 
g/K can then be obtained by making appropriate choices 
for (7c, ri, and Ai. The only conditions that we need 
to satisfy are Ai ^7,- and gc,^ ^ Ai. That allows us 
some flexibility to pick a smaller gc, which determines the 
cavity length, and a larger fl, than is the case for the adi- 
abatic approach proposed in [s^. For instance, using the 
hyperbolic secant pulse shape and the optimum control 
pulse shape, given by Eq. (^5)1 . we have gc^{t) = g'^Q'{t), 
where n'{t) is only a function of kT and not gc- There- 
fore, for a fixed value of kT, the only way to increase 
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FIG. 7: (Color online) The loading probability for the adi- 
abatic method as a function of hiT versus gc/{KAi), at the 
two-photon resonance Ai — A2 with (j}z{t) = 0, 7r = 0, and 
several values of kT. These curves were obtained by numeri- 
cally solving Eq. (|18|) using the control pulse shapes, shown in 
the inset, that correspond to a hyperbolic secant input pulse 
shape. 



the coupling rate in the adiabatic case is to use a shorter 
cavity, which yields a larger gc. 

For nonzero values of 7r/Ai, there is another param- 
eter that we can vary to optimize the loading proba- 
bility in our non-adiabatic scheme: it is gd^, which 
appears in the exponents —[g'^jr/^i + K±]{t — t) = 
-[K + g{gc/n + n/gc){ir/Ai)±^]{t-T)/2mEq. To 
increase the loading probability, for fixed values oi g/n 
and 7r/Ai , we need low values of these exponents. These 
can be achieved at gd^ ~ 1- The exact optimum value 
for gd^ that minimizes the decay effect also depends on 
the values of 7 and ^. For the sake of comparison, how- 
ever, we assume that 7r/Ai is sufficiently small that we 
can neglect these subtleties, and we only focus on the 
interplay between the coupling rate and the input band- 
width. So, for the rest of this section, we assume 7^ = 0. 

It is interesting to find the dependence of the loading 
probability, for the adiabatic scheme proposed in [33| . 
on the coupling rate gc- For this purpose, we have nu- 
merically solved the Schrodinger equations in Eq. (jlSp . 
for r2(i) given by Eq. at the two-photon resonance 

Ai = A2 and for <f)z{t) = 0. In Fig. [7] we have plotted 
|ce(5T)p versus an effective coupling rate g' = g^/Ai. 
Choosing t = 5T ensures that the loading process has 
ended, and therefore |ce(5T)p is effectively the loading 
probability. This figure shows that for the two-photon 
resonance case, higher gc values yield higher loading 
probabilities. The effective coupling rate g' that we need 
is about 2k for 90% loading probability at kT = 5. It 
can be seen that there is an advantage to using pulses 
with higher values of kT, because our adiabatic scheme 
is more efficient for longer input pulses. 

We can also solve the Schrodinger equations in 
Eq. m]), for n{t) given by Eq. at zero effective 

detuning, i.e., when (p^it) = Ai - A2 [g^ - f^^(t)]/Ai, 
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FIG. 8; (Color online) The optimum coupling rate and 
the maximum loading probability achieved at this rate for 
both adiabatic and non-adiabatic approaches. (Non-adiabatic 
curves start from kT — 0.5; adiabatic curves start from 
fcT = 4.5.) For all curves, we assume that A = 7^ = 0. A 
hyperbolic secant pulse shape has been used for the incoming 
photon. 



as proposed in [2^. Note that the control field's Rabi 
frequency given by Eq. (j23[) has been obtained under the 
two-photon resonance and dark-state conditions, and it 
is therefore not the optimum pulse shape for the zero- 
effective-detuning case. It can be verified, numerically, 
that this control field provides a close-to-optimum per- 
formance. This is in accord with the results reported in 
[3^ implying that the loading probability will only be 
slightly affected by small deviations from the optimum 
control field pulse shape. For zero-effective detuning, in 
analogy to the non-adiabatic case, we observe the exis- 
tence of an optimum value for g' /k. Figure [8] shows the 
optimum coupling rates and the maximum loading prob- 
abilities, Pioadingi achieved at these rates as functions of 
kT for the adiabatic and non-adiabatic approaches us- 
ing a hyperbolic secant pulse shape. Here, {g' / Kjopt is a 
near-unity constant, for all values of kT, whereas {g/ n)opt 
decreases toward zero with increasing kT . The maximum 
loading probability that can be achieved by the adiabatic 
approach asymptotically goes to one as kT increases. Its 
value is about 75% for kT = 4.5, but we should bear in 
mind that the system will be more vulnerable to sponta- 
neous decay when we use long pulses. 

Because the non-adiabatic approach has no constraint 
on input pulse width, it allows a larger class of single- 
photon sources to be used in quantum communication 
systems, which can result in higher throughputs. For in- 
stance, consider system which uses a downconverter 
single-photon source at kT — 1, and non-adiabatic load- 
ing at the optimum coupling rate. We will compare it 
with system B, which uses adiabatic loading with the 
same source but at kT — 4, the minimum required for 
the adiabatic case in our examples. Both systems use 
identical quantum memory modules, so their k values 
are the same. The total throughput is proportional to 
the source rate and the loading probability. The source 
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FIG. 9: (Color online) Loading probability versus normal- 
ized timing offset in the photon arrival. In the adiabatic case 
(the left curve), this timing error results in applying the con- 
trol field earlier or later than the correct time. In the non- 
adiabatic case (the right curve), it results in stopping the 
control field earlier or later than T^oad- Both curves assume 
A = 7r = for their respective optimum coupling rates for 
hyperbolic secant input pulse shapes. 

rate for system A is four times that of system B, be- 
cause we are allowing a larger output bandwidth in the 
former case. From Fig. [51 the maximum achievable load- 
ing probability for system A is above 75%, while that for 
system B is at most one. It folUows that system A can 
have a 4 X 0.75 = 3 times higher througput than system 
B. 

Finally, Fig. [9] shows the sensitivity of both schemes 
to timing offsets. Here, Toffsct refers to the time that we 
stop the control field in the non-adiabatic approach, and 
it represents a time shift in the adiabatic approach, i.e., 
using f2(i — Toffsct) instead of fl(t). Both cases may occur 
if we have an inaccurate estimate of the photon arrival 
time. This figure shows that both schemes have almost 
the same tolerance for timing offsets. A 50% loading 
probablity (3 dB loss) is achievable even if we are about 
±T/2 off from the correct loading time. Although it is 
not shown in this figure, the same result holds if we use 
the adiabatic approach under the two-photon-resonance 
condition. 



C. Double-A Trapped Atoms 

The results we obtained above are directly applicable 
to several other problems. For instance, consider the 
double- A atom shown in Fig. llOf a). Here, we are trying 
to load the memory with a single photon in an arbitrary 
polarization state: 

li^o) - J du;4>{u;)ia\l^),^+l3\U,_), (26) 

where jap + — 1, and |laj)cr± refers to a single- 
photon state at frequency lo in the <7± (right/left cir- 
cular) polarization. The goal of the loading process is 
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FIG. 10: (Color online) (a) Loading a double-A atom driven 
by a polarization-state single photon. By adiabatically elim- 
inating the upper states, this system becomes equivalent to 
the V^- level system shown in (b). (c) Reduction of the level 
loading problem to two two-level loading problems by means 
of linear superposition. 

to absorb this photon, by properly choosing z-polarized 
control fields — according to either the adiabatic or non- 
adiabatic mechanisms — and store its polarization infor- 
mation in the metastable levels je-i.) and |e_) that corre- 
spond to the (T_|_ and cr_ polarizations, respectively. The 
loading probability is therefore determined by the prob- 
ability of ending up in |e) = a\e^) + /5|e_). 

In order to calculate the loading probability, we can use 
the results of the previous subsection to reduce each A- 
configuration leg of the atom to an equivalent two-level 
system, as shown in Fig. fTOF b). The resulting config- 
uration is a I^-level atom, illuminated by an arbitrary- 
polarized single photon. This problem can be solved 
by noting that the full Schrodinger equations — for the 
state of the reservoirs, cavity modes, and the atom — 
that model the system evolution are linear in response to 
the superposition state in Eq. ([26]). Therefore, each leg 
of the I^-level atom is only driven by the correspond- 
ing part in the initial state, as shown in Fig. Iior c). 
In fact, for a / c?w0(ct;)|lij)o-+ as the input, C-{t) — the 
probability amplitude of being in the state |e_) — will 
be zero. We only get a nonzero c_(i) from the input 
(3 J duj(j){uj)\luj)a-^, in which case the C-{t) evolution co- 
incides with that of a two-level system. For instance, if 
we are using the non-adiabatic loading mechanism, c_ (t) 
is given by Pce{t), using Ce{t) from Eq. (^5)) . The same 
argument holds for c+(i), so that the overall probability 
of ending up in |e) is (jap |/3|2)|ce(i)P = |ce(i)p. This 



proves the equivalence of loading a y-level atom to the 
loading of a two-level atom when both are driven by sin- 
gle photons. This result also extends to cases in which the 
two legs of the atomic system are nonlocal — e.g., when 
a pair of two-level trapped-atom memories are driven by 
a photon-number entangled state. Such a system is an 
example of the to-the-memory configuration for entan- 
glement distribution, and it can be used in teleportation 
systems. 



III. MIT-NU ENTANGLEMENT 
DISTRIBUTION 

With the results that we have developed for trapped- 
atom quantum memories in the previous section, we can 
address the performance of quantum communication sys- 
tems that use this type of storage. Our focus on this 
section will be on the MIT-NU architecture for quantum 
communication [2^ . The MIT-NU construct uses the to- 
the-memory configuration in Fig. [IJa) to distribute and 
share polarization entanglement between two rubidium 
atoms. Because the rubidium atoms can be treated as 
double- A quantum memories [11], MIT-NU entanglement 
distribution is essentially a loading problem, in which we 
are interested in transferring an entangled state to two 
trapped atoms. All previous analyses of the MIT-NU 
architecture [lES, EH, however, have employed a cold- 
cavity approach in which each optical cavity — that would 
hold an *^Rb atom in the actual implementation — is re- 
garded as empty, and the loading probability is calculated 
by determining the probability that the state of the in- 
tracavity photon fields at the end of a loading interval is 
the desired singlet. In this section, we provide the first 
loading analysis of the MIT-NU system that includes the 
presence of atoms within the quantum memory units. 

FigurefTlTa) shows a schematic of the MIT-NU system: 
QMi and QM2 are trapped rubidium atom quantum 
memories, each Lq km away — in opposite directions — 
from a dual optical parametric amplifier (OPA) source. 
Each optical parametric amplifier in the dual-OPA source 
is a continuous-wave, type-II phase matched, doubly- 
resonant amplifier operating at frequency degeneracy 
Us — to I- Its signal (S) and idler (J) outputs comprise a 
stream of orthogonally-polarized photon pairs that are 
in a joint Gaussian state [23|. By coherently pump- 
ing two of these OPAs, 7r-rad out of phase, and com- 
bining their outputs on a polarizing beam splitter, as 
shown in Fig. Illf b). we obtain signal and idler beams 
that are polarization entangled [Tj. Moreover, double- 
resonant operation greatly enhances the brightness of 
these beams, making this entanglement source compati- 
ble with loading atomic memories [6], something that is 
not the case for standard ~THz bandwidth parametric 
down-converter sources. The signal and idler beams are 
routed down separate optical fibers to the trapped-atom 
quantum memories. To do so efficiently, their wavelength 
is chosen to lie in the low-loss 1.55 /xm band. However, 
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FIG. 11: (a) MIT-NU architecture for long-distance quan- 
tum communications consisting of a dual-OPA source that 
produces polarization-entangled photons, and two quantum 
memories, QMi and QM2, separated by 2Lo km. (b) Dual- 
OPA source of polarization-entangled photons. OPAs 1 and 2 
are coherently-pumped (yr-rad out of phase), continuous- wave, 
type-II phase matched, doubly-resonant amplifiers operated 
at frequency degeneracy whose orthogonally-polarized signal 
({Sfc}) and idler {{Ik}) outputs are combined, as shown, on 
the polarizing beam splitter (PBS), (c) Notional schematic for 
the relevant hyperfine levels of *^Rb. Each quantum memory 
consists of a single trapped rubidium atom that can absorb 
arbitrarily polarized photons, storing their coherence in the 
long-lived D levels. A non-destructive load verification is ef- 
fected by means of the A-to-C cycling transition. 



assume that both memories have stored photon coher- 
ences and so are ready for the rest of the teleportation 
protocol, i.e.. Bell-state measurement, classical commu- 
nication of the result, and single-qubit rotation [TT'|. No 
entanglement purification step was envisaged in the orig- 
inal description of the MIT/NU architecture, although 
one certainly could be added if multiple pairs are stored 
at each location. 

In order to find analytical results for the MIT-NU load- 
ing problem, we make several simplifying assumptions. 
First, we replace the OPA devices in the entanglement 
source with simple spontaneous parametric downconvert- 
ers (SPDCs). This is actually how most realizations of 
this system work [2^ Is^] , and, except for their flux and 
their output bandwidths, SPDC systems have the same 
physical characteristics as OPA sources. The second as- 
sumption is the biphoton approximation to the SPDC 
output. The latter is mostly a vacuum state, plus a 
small biphoton component, and a much smaller multi- 
pair contribution [23|. The vacuum term can be eas- 
ily recognized by our non-destructive loading verification 
technique, and can therefore be ignored. The multi-pair 
case is of minor concern in the context of loading, be- 
cause its degrading effect has been accounted for in the 
previous analyses [3J|. So, in this section, we only con- 
sider the case in which the type-II phase-matched SPDC 
output, operating at frequency degeneracy, is a biphoton 
state in the following general form [43| 



dujs / (itJ/(/)(w5,cj/)|l 



(27) 



because the quantum memory makes use of the *^Rb line 
at 795 nm, quantum-state frequency conversion is used, 
at each memory location, to convert the polarization en- 
tanglement from 1.55 /im to 795 nm [42]. As explained in 
[13], additional steps are taken to, prior to the upconver- 
sion, to compensate for any polarization transformation 
encountered in transmission from the source to the mem- 
ory. This is accomplished by monitoring the polarization 
change incurred on a strong pilot pulse and using that 
information to drive polarization controllers that com- 
pensate the signal and idler beams. 

A schematic of the relevant hyperfine levels of ®^Rb 
is shown in Fig. IllT c). The memory atoms are initially 
in the ground state A. From this state they can absorb 
a photon in an arbitrary polarization and, by means of 
a Raman transition, transfer that photon's coherence to 
the long-lived D levels for subsequent use. However, be- 
cause propagation and fixed losses may destroy photons 
before they can be stored, and because both memories 
must be loaded with a singlet state prior to performing 
qubit teleportation, the MIT-NU architecture employs a 
clocked loading protocol in which the absence of fluores- 
cence on the ^-to-C cycling transition provides a non- 
destructive indication that a memory atom has absorbed 
a photon. If no fluorescence is seen from either the QMi 
or QM2 atoms in a particular loading interval, then we 



where J dujs J dujj\(l)(ujs,uji)\'^ — 1, and o and o represent 
two orthogonal polarizations. For the particular config- 
uration in Fig. [TlT b). the desired output state will then 
be 



li'out) = / dujs / duji4>{ujs,i^i)\-ip^s,'^i) 



SI' 



(28) 



where 



/S.cr 



SI 



V2 



(29) 

is the singlet state, which is invariant with respect to 
polarization basis. Finally, we assume that there is no 
loss and no delay in the channel so that the preceding 
state is the input state to the quantum memories. 

Figure [T2l^a) shows a schematic of the MIT-NU load- 
ing problem with the above considerations taken into ac- 
count. Here, similar to the previous section, we assume 
that the signal and idler photons are spatially matched 
to one of the cavity's spatial modes with center frequency 
luq. Depending on their polarizations, they can drive one 
of the transitions \g) — *■ |s) (for a+-polarization) and 
Iff) ^ k) (fo^' c_ -polarization) with time- independent 
coupling rate gc and detuning Ai. The second transition 
is facilitated by a 2-polarized classical field, which induces 
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FIG. 12: (Color online) (a) Loading a pair of double-A atoms 
illuminated by a polarization-entangled biphoton. This is 
a fair approximation to the MIT-NU loading problem, (b) 
Breaking the loading problem in (a) into two simpler loading 
problems for a pair of A-level atoms, each illuminated by a 
biphoton state. 



a Rabi frequency with detuning A2. The goal of mem- 
ory loading is to absorb the entangled photons and store 
their coherence in the metastable levels of each memory, 
i.e., to end up in the state |e+)s|e_)/ — |e_)s|e+)7, where 
the subscript S/ 1 refers to the atom that is driven by the 
signal/idler photon. 

The above problem can be reduced to the two simpler 
problems, shown in Fig. I12f b). by using the linearity in 
the input superposition state. This is the same technique 
that we used in Section HTC] to reduce a T^-level atom to a 
pair of two-level atoms. Now that our problem has been 
reduced to the loading problem for two A-level atoms, 
we can employ the adiabatic and non-adiabatic loading 
mechanisms described in Section fll Bl We provide ana- 
lytical results for the loading probability associated with 
either of the systems shown in Fig. [TWb) , which can be 
directly applied to the case shown in Fig.fT^a). or equiv- 
alently, to our approximation to the MIT-NU loading 
problem. Because the two subsystems in Fig. [TUJb) have 
the same loading behavior — independent of the polariza- 
tion of the incoming photons — we simplify our notation 
by omitting the polarization information. The driving 
state will then be written as follows 

IV'o) = y" dujs J dLJi(j){uJs,^^i)\'^i^s)s\i-i^i)i: (30) 

where J dujs J dujj\(l)(ujs^^i)\'^ — li ^nd the subscripts S 
and / refer to the two independent reservoirs that inter- 
act with our memory cells. 

In this section, we only derive the loading probability 
for our non-adiabatic scheme. As we showed in Eq. 
if the decay rates from the upper states are much smaller 
than the employed detuning, each A-level atom can be 
approximated by a two-level atom with an effective cou- 
pling rate g = gci^/Ai and a detuning A that can be 




FIG. 13: (Color online) A two-level-atom model for the sys- 
tems shown in Fig. I12f b'). Here, we have assumed that the 
Rabi frequency fi is a constant and we have adiabatically elim- 
inated the upper state in the A-level atoms. In the absence of 
spontaneous decay, the effective coupling rate is g = gcf^/Ai. 
The effective detuning A can be made zero by a proper choice 
of parameters, as shown in Eq. (fT8|l . 



made zero by a proper choice of parameters. This re- 
duces the loading problem in Fig. flW b) to the one shown 
in Fig. 1131 in which a pair of trapped two-level atoms 
are illuminated by a biphoton state j-^o) as given by 
Eq. ([50)1 . Here, the loading event corresponds to pop- 
ulating both atomic excited states, i.e., to ending up in 
the state |e)s|e)/. 

To analyze this loading problem, we can again use the 
superposition trick from the previous subsection to find 
Cee(i), the slowly- varying probability amplitude at time t 
associated with |e)s|e)/. Using 0^1^,10 At) to denote the 
slowly-varying probability amplitude at time t associ- 
ated with |e)5|e)/ in response to the two-photon state 
0(t^s,w/)|l<^s)s|lwj)/, we have 



Ceeit) = I duJs I duJlCu,s,uj,it), (31) 

where, from Eq. (jlip . 



c^s,ujiit)dujsduji 



2 pt pt 



Kg 

^ g-i(a;s-wo)Tg-i(a)7-a;o)-r' 



X (e 
X (e 



dr I dr' du)sdiOi(j){u)ST'-^i) 



)■ (32) 



-K+(t-T') _ -K_(t-r') 



The above equation results in the following expression 
for the time-dependent probability amplitude for being 
in the state |e)s|e)/: 



Cee{t) 



^ f dT f dr'Mry) 



X (e' 



-K+(t--T') _ -fc_(t-r') 



), (33) 



where 



$fc(r,T')= jdus ydu;/e-^("^-"°)--^(-^-"'')^>(c^s,^/). 

(34) 

Note that in the above integral only the symmetric part 
of the pulse shape, [$h(r, r') -I- <i>f,(r', r)]/2, results in 
a nonzero value for Ceeit). Also, it has been implicitly 
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assumed that the pulse shape <I>6(r, r') is nonzero only 
for T, t' > 0. Otherwise we have to change the lower 
limits of the above double integral. 

A. Numerical results 

Here, we present some numerical results for the per- 
formance of MIT-NU non-adiabatic loading mechanism. 
Our goal is to find the dependence of the loading proba- 
bility on the bandwidth of the input pulse as well as on 
the atom-light coupling rate. This loading probability is 
given by |cee(i)P ^ obtained in Eq. ([55]) . The input pulse 
shape that we consider here corresponds to the output of 
an SPDC, which is given by [i^l 

4}{ujs,uji) = A4>+{ujs,LJi)(l)^{ujs,Lui). (35) 

In this expression, ^ is a normalization factor, and 



ns{uJs)ni{uJi) 
ns(wo)n/(wo) 



(/>p(ws+w/), (36) 



where jt-s// is the downconversion crystal's refractive in- 
dex for the signal/idler beam and (f>p{uj) is the pump 
spectral pulse shape. The other term in Eq. ([55]) is the 
phase-matching function. 



sin[Afc(ws,w/)i/2] 



(37) 



where L is the crystal length, and A/c(w5, ujj) = kpiujs + 
^i) ~ ks(j-^s) — ki{uji), with kg/j/p being the wave 
number of the signal/idler/pump beam. We assume 
that the crystal is phase matched at degeneracy, i.e., 
kp{2u!o) = ks{uJo) + fc/(wo), and that it also satisfies the 
extended phase-matching condition, k'p(2LUo) — A;^(wo) = 
fcj(cjo) — k'p{2iLjo) = Sk, where k'j({uj), X — S,I, P, is the 
derivative of kx with respect to w. A first-order linear 
approximation will then yield Afc(u;5,w/) « {los — (^i)Sk. 
Consequently, (p±{u}s,'^i) is only a function of ± w/. 
The temporal baseband pulse shape needed for our load- 
ing probability calculation is therefore 



A' 



t + u 



t — u 



(38) 



where A' = ylwo/[«s(t^o)'^/(^^o)], ^b{t) = 
e^*"^"* / dw(/)p(a;)e~''^* is the baseband pump pulse 
shape, and $_(i) = {2'K/5k)[u{t + To/2) - u{t ~ To/2)], 
where u{t) is the step function and Tq — 5k L. 

Figure [Ml shows the loading probability for the archi- 
tecture in Fig. [13] when hyperbolic secant pulse shape 
^b{t) = y27rsech[4(t - To - 2T)/T] is employed. We 
have shifted the pump pulse shape by 2T-\-Tq so that, ef- 
fectively, ^h{t^ u) = for t < or M < 0. Now, in addition 
to its dependence on g/n and kT, the loading probabil- 
ity also depends on kTq, which represents the bandwidth 



e 
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1-^ 
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FIG. 14: (Color online) (a) MIT-NU loading probability for 
different values of g/n at nT = kTq — 2. The maximum 
probability in each case is a function oi g/n, which implies the 
existence of an optimum value for g/n. (b) MIT-NU loading 
probability for different values of kTq &t g/n = 1 and kT = 2. 
In all curves, we have used a hyperbolic secant pulse shape 
for the pump beam with 7' = 0. 



of our downconversion process. In Fig. [Hf a). we have 
plotted the loading probability for different values of g/ k 
with fixed values of kT and kTq at 7' = 0. Similar to 
what we found in the single-atom case, this figure shows 
that there exists an optimum value of g/ k whose corre- 
sponding peak loading probability is maximum. Hence, 
if we turn the control field off at this particular time, the 
chance of finding both atoms in their excited states is 
maximum. Figure I14f b) shows the same property from 
a different perspective. It plots the loading probability 
with kT and g/ k fixed for several values of kTq. Here, 
we see that there is an optimum value of kTq for which 
the loading probability is maximum. 

Figure [T51 shows the optimum value oi g/ n and its cor- 
responding loading probability as a function of kT and 
kTo at 7' = 0. From Fig. [TSTa'). we observe that higher 
values of kT or kTq, lead to lower values of g/n. This 
is another manifestation of the trade-off between the in- 
put pulse bandwidth and the required coupling rate. The 
lower the bandwidth is the lower the optimum coupling 
rate will be. Note that this regime of operation is not 
applicable to the adiabatic-passage approach on the two- 
photon resonance, in which the higher the coupling rate 
is the higher the loading probability will be; see Fig. [T] 
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IV. CONCLUSIONS AND DISCUSSION 




FIG. 15: (Color online) (a) The optimum value oi g/n ver- 
sus kT and kTo for the MIT-NU non-adiabatic loading prob- 
lem, (b) The corresponding maximum loading probability at 
{g/fi)opt- We have used a hyperbolic secant pulse shape for 
the pump beam with 7' = 0. 



As a result, the non-adiabatic loading mechanism has 
an advantage in that it allows us to use larger cavities 
and higher input bandwidths. The price, however, is the 
value of loading probability that can be achieved. As 
shown in Fig. I15f b). for non-adiabatic loading this prob- 
ability is typically on the order of 70%-80% — rather than 
100%. From Fig. [T5f b). we see that higher loading proba- 
bilities occur for kTq, kT > 2, where the cavity linewidth 
is wider than the bandwidths associated with the input 
pulse. 

The initial plan for the MIT-NU architecture called for 
cavities with k « 5 MHz and gc ~ 10 MHz. That can be 
achieved with a cavity of length ~ 400 fim [l3] , which 
corresponds to a Purcell factor Fp = 3c/(47r^KLc) ~ 
4 X 10*. That leaves a considerable flexibility for choos- 
ing n and Ai to operate at the optimum coupling rate 
g = gc^/ Ai. On the other hand, in order to have 
a high-throughput system it is required to have opti- 
cal sources with narrow bandwidths comparable to k. 
This can be achieved with the recent progress in build- 
ing cavity-enhanced SPDCs @. Overall, although the 
loading probabilities that we have obtained in this section 
are under somewhat idealized conditions, their values are 
high enough to make the 5dB fixed loss used in [23, 3^, 
not an unduly optimistic, but a conceivably conservative 
assumption. 



Many proposals for entanglement distribution over 
long distances rely on our ability to transfer the state of 
a single photon to an atomic quantum memory. Such a 
memory unit is comprised of either a single atom trapped 
in an optical cavity or an ensemble of atoms. In this 
paper, we addressed the problem of loading a trapped 
atom with two-level, A-level, and double-A-level config- 
urations. The general approach suggested for solving 
this loading problem is based on adiabatic-transfer tech- 
niques. In this paper, we looked at the problem from 
a different viewpoint and realized that a non-adiabatic 
mechanism may have characteristics that are compara- 
ble to, and in some aspects even better than, those of 
adiabatic loading. Our non-adiabatic method is based 
on the Rabi oscillation that naturally occurs when a sin- 
gle photon illuminates a trapped-atom quantum memory. 
We then freeze this oscillation when the probability of be- 
ing in the desired state is maximum. We showed that by 
properly choosing the system parameters, e.g. atom-light 
coupling rate versus input bandwidth, we could achieve 
loading probabilities above 90% in memory units with 
A or double- A atomic configurations. In terms of tim- 
ing sensitivity, we showed that both adiabatic and non- 
adiabactic mechanisms have the same tolerance to errors 
in their estimate of photon's arrival time. 

Whereas the adiabatic passage is inherently suited for 
long input pulses, our non-adiabatic technique does not 
impose any restrictions on the source bandwidth and 
therefore allows a larger class of optical sources to be 
employed. That freedom can potentially increase the to- 
tal system throughput for distributing entanglement. 

For any input photon bandwidth, we realized the exis- 
tence of an optimum coupling rate, in our scheme, that 
maximized the loading probability. This optimum cou- 
pling rate slightly increased when we introduced spon- 
taneous decay, but it nevertheless remained in a range 
achievable by current technology. We also identified 
a trade-off between the required atom-light coupling 
strength and the input bandwidth, which showed that 
in order to accommodate a broadband source we needed 
to enhance our coupling rate. Conversely, the optimum 
coupling rate approached zero, although very slowly, as 
we decreased the bandwidth associated with the driving 
photon. We observed the existence of an optimum cou- 
pling rate in some adiabatic loading techniques as well. 
They did not, however, offer such flexibility in working 
either with lower coupling rates or shorter input pulses. 

Finally, we applied our initial results to the MIT-NU 
loading problem. The MIT-NU architecture uses a pair 
of trapped- rubidium- atom quantum memories, which are 
illuminated by the signal and idler outputs from an ultra- 
bright doubly-resonant dual-OPA source of polarization- 
entangled photons. By approximating these outputs as 
a general biphoton state, we derived, analytically, the 
loading probabilities for the non-adiabatic approach. In 
this case, we again observed the existence of an optimum 
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coupling rate as a function of the driving pulse param- 
eters. We showed that loading probabilities above 80% 
were achieveable at these optimum coupling rates, pro- 
vided that the input pulse's bandwidth is narrower than 
that of optical cavities. 

Although all our calculations are based on trapped- 
atom quantum memories, most of the results presented 
here can be extended to cold atomic-ensemble memories 
by including an additional multiplicative factor \fN , with 
N being the number of atoms in the ensemble, in the cou- 
pling rate g. This factor accounts for the collective effect 



of ensemble. For room-temperature atomic ensembles, 
one should also account for the inhomogcneous broaden- 
ing due to atomic motion, whose treatment needs a more 
general model than what we described in this paper. 
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